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Abstract 

We consider a cosmological model with a variable gravitational constant, G, 
based on a scalar-tensor theory. Using the recent observational data for the Hub¬ 
ble diagram of type la supernovae (SNela) we find a phenomenological expression 
describing the variation of G. The corresponding variation of the fine structure 
constant a within multidimensional theories is also computed and is shown not to 
support known constraints on Aa/a. 


1 Introduction 

The analysis of the observations of the Hubble diagram of distant type la supernovae 
Pi21iH| provide evidence that the universe has been accelerating recently, at z < 0.5, 
and decelerating at earlier stages Pi- The Friedmann cosmology without cosmological 
constant and with zero curvature — as indicated by the recent CMB Boomerang and 
Maxima experiments Pi — cannot explain such an evolution of the universe jZj. The 
accelerated behaviour can be attributed to a “dark energy” with negative pressure, the 
simplest possibility being the introduction of the cosmological constant in accordance with 
cosmic concordance model Qm ~ 0.3, Da ~ 0.7 — see jS] for a review. 

An alternative solution is to modify the gravitational theory, for example, by allowing 
a time variation of Newton’s constant. The possibility of a time variation of fundamental 
constants of nature, in particular of the fine structure constant, a, and of the gravitational 
constant, G — first considered by Dirac in the framework of his Large Number hypothesis 
[H] and later developed in ^ within an alternative theory of gravitation (see references 
CH and [12] for more details) — has been recently a subject of numerous studies (see 
for example references [m ca d [13 for recent reviews and extensive bibliography). 
It is worth mentioning that many theoretical approaches, such as models with extra 
dimensions, string theories or scalar-tensor models of quintessence, contain a built-in 
mechanism for a possible time variation of the couplings. Astronomical measurements 
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allow to constrain such hypothetical variations. As a matter of fact, local constraints on 
the rate of variation of G can be derived, for example, from Lunar Laser ranging cniini, 
whereas constraints at cosmological distances can be derived, amongst other methods, 
from the Hubble diagram of distant SNela [TH] , 

On the other hand, in the framework of models with a varying gravitational constant 
it would be valuable to get a phenomenological expression for its variation. That is, 
to obtain an approximate form of the function G{z), where z is the redshift. Getting 
such description is the goal of the present paper. We would like to note that htting 
the Hubble diagram of SNela within models with a variable gravitational constant for a 
particular parametrization and A = 0 was studied in ng. Using the considerably better 
observational data available now j3] which extend to much larger distances we will hnd 
a more accurate approximation for G{z) for a larger interval of look-back times. It is 
worth mentioning as well that a procedure of reconstruction of a general scalar-tensor 
model (the scalar held potential and the functional form of the scalar-gravity coupling) of 
dark energy from cosmological observational data, particularly the luminosity distance, 
was hrst developed in ra. 

Once the phenomenological form of G{z) is obtained it can be compared with pre¬ 
dictions of cosmological models and/or contrasted with other astrophysical observational 
constraints. For instance, models with extra dimensions incorporate a natural mechanism 
for the space and time variation of various fundamental constants, which was appar¬ 
ently studied for the hrst time in ^I] and later on in a number of articles. It should be 
also mentioned that in j22j the relation between the time variation of the hne structure 
constant and that of the gravitational constant was studied for three classes of theories 
— namely, for the pure Kaluza-Klein theory, for Einstein-Yang-Mills theories and for 
Randall-Sundrum type models. Using the relation between a and G in a given model 
and the phenomenological expression for G{z) one can then obtain an estimate of the 
variation of the hne structure constant. This prediction can therefore be then contrasted 
with the observational constraints on the variation of a, which has recently been a subject 
of intensive studies. Particularly, using the many multiplet method it has been claimed 
1231211123 that the hne structure constant a was smaller in the past. However, a similar 
analysis carried out in [23 and in m using a diherent hne htting code and data sample 
of better quality shows that the measurements are consistent with zero variation within 
the observational uncertainties and, consequently, these results do not support the claims 
by previous authors. 

The plan of the paper is the following. In Sect. |21we outline a theoretical scheme for 
a variable gravitational “constant” and derive a generalization of the Hubble law for this 
case. A phenomenological description of the function G{z) which hts the Hubble diagram 
of SNela is then found in Sect. El In the next section the correlated variations of the hne 
structure constant and of the gravitational constant for models with extra dimensions 
is discussed. Finally, in Sect. El our main conclusions are presented, followed by some 
discussion of our most important results. 
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2 Variation of G in scalar-tensor theories 


Theories of gravity in which the gravitational “constant” G varies with time and cos¬ 
mological models based on them have been extensively studied in the literature (see, for 
example, j2H] and references therein). One of the most natural and relativistic covariant 
ways to describe the variation of the gravitational constant is to interpret it as a scalar 
field 0. This can be done self-consistently in the framework of scalar-tensor theories of 
gravity of the Jordan-Brans-Dicke type Iinii2i with the action given by 

S = J + IGTrL^^ , (1) 

where the function w = w{(j)) determines the coupling between the scalar field and gravity. 

In considering the cosmic evolution of the scale factor a{t) of the Friedmann-Robertson- 
Walker metric and of the scalar field 0 in Eq. ([Q) we assume for simplicity that w is con¬ 
stant. Then, the Hubble parameter H = a/a is determined by the Friedmann equation 

I2BI 




/a\^ Svr k 0 a tc 0^ A 

Va/ 30^ a^ 0 a 6 0^ 3 ’ 


( 2 ) 


where A is the cosmological constant and k is the curvature parameter. We furthermore 
assume that the universe contains a simple perfect fluid described by the equation of state 


P = (7 - 1)P- (3) 

Eqs. m and © and the energy conservation condition p + S'jHp = 0 have to be comple¬ 
mented with the acceleration equation for 0 j28j : 

In what follows we will consider a and 0 to be functions of the redshift 2 ;. To convert time 
derivatives to the derivatives with respect to 2 ; we use the standard relation: 


A 

dt 


dz 


Denoting the 2 ;-derivatives with prime we get relations of the type 0 = —H{1 -|- z)(l)'. 

By considering the weak-field limit in the scalar-tensor theories the following relation 
between the gravitational constant and the scalar held 0 can be established UDI 


G(z) 


4 -I- 2ta 1 
3 -|- 2ta 0 ( 2 ;) 


(5) 


Using these expressions, the Hubble law — given by Eq. (j2I) — can be written after 
some algebra in the following form 






9{z) 


( 6 ) 
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The function g{z) is given by 

g{z) — 1 + {1 + z)—-—{1 + z) ! (^) 

where Go = G'(O) is the present-day value of the gravitational constant, go = g{0), and 
the parameters CIm, and are related to the standard ratios 


SttGoPo 
through the following relations 

3 T ‘2w 


ao-no 


r^A = 


From Eq. 


— 

it follows that 


, = ^, Q 

go 4 +2m; go 


A — 


A 


9o ' 


( 8 ) 


+ ^A — 1- (9) 

We would like to note at this point that a particular case of Eqs. was discussed 

in pnj . 

Finally, the luminosity distance is calculated via the standard formula, which in 
the flat case has the form 


rz flqi 

dL = c{l + z) (10) 

Jo H [u) 

Similar to the way in which it was done in , for the calculation of the Hubble diagram 
of SNela we will use the Chandrasekhar mass model for the SNela light curve. According 
to this, the peak luminosities of SNela are proportional to the Chandrasekhar mass (L oc 
Mch) and, therefore, scale as L oc This result has been validated by detailed 

numerical calculations of exploding white dwarfs. As a result the apparent magnitude is 
given by 


/ \ I 1 .J- n dL{z)Ho __ 15 G{z) 

m{z) = Mo+ 5 log-^ + 25 + — log —— 

c 4 Go 


( 11 ) 


where Mq is the absolute magnitude. We will use Eqs. (fTn|l and (FTTll to fit the Hubble 
diagram of distant SNela with a certain parametric representation for G{z). 


3 A phenomenological description of the variation of 
the gravitational constant 

As it has been already commented before, our goal is to obtain an empirical description 
of the variation of the gravitational constant as inferred from the observational data of 
SNela, including the most recent and reliable datasets j2]. For this purpose we use a 
simple phenomenological expansion of the function G{z) in powers of z: 
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G{z) — Go(l + Pi^ + P 2 z‘^ + Psz^ + (12) 

Such approximation is in the spirit of phenomenological descriptions of the scale factor, 
deceleration parameter or equation of state emloyed previously by other authors pin ED- 
The hrst coefficient in Eq. m, Pi, is determined from experimental bounds on the 
time derivative of the gravitational constant at the present time {G/G)o = {G/G)t=t^^^- 
For convenience we translate these bounds in the ones on r; = G'{Q)/Gq using the relation 


g'(0) 

Go 




For this estimate we will take the present value of the Hubble constant to be Hq = 
63 km/s/Mpc. 

There are a number of constraints on {G/G)q obtained from very different observations 
and methods — see, for instance, m and na and references therein. For example, 
the Lunar Laser ranging experiments yield |G/G|o < 8 x 10“^^ yr“^ j1 61 IT7j . whereas 
the improved constraints on the post-Newtonian parameters give an upper bound of 
10“^^ yr“^ [22] • Summarizing, one can see that the parameter p can take values satisfying 
roughly \p\ < 0.01, which is similar to the estimate obtained in Actually, as far as 
p is small enough, the values of the other coefficients depend very weakly on its precise 
value. Moreover, we have checked that our hnal result is not sensitive to the value of 
p within the interval —0.01 < p < 0.01. Since p = pi the bound on the local rate of 
variation of G determines the linear term in Eq. m- For the forthcoming analysis we 
adopt Pi = p = —0.01. 

In what follows we truncate the expansion given in Eq. m at a certain order, substi¬ 
tute this polynomial expression into Eq. m, calculate the apparent magnitudes of SNela 
as given by Eq. m in terms of the coefficients pi and compare them with the Hubble 
diagram based on the observational data of SNela. 

Let us consider the case of the flat and matter dominated universe without the cos¬ 
mological constant. Namely, we set = Ha = 0 and 7 = 1 in the Hubble law, Eq. 

An important observation is that with such assumptions the value of the ratio Qm turns 
out to be hxed by the value of pi. Indeed, from m one gets 


A + 2w A + 2w 
~ 3 + 2w^° ~ 3 + 2w 


(i+Pi-|ri). 


(13) 


The value of the Brans-Dicke parameter w depends on the specihc model. For exam¬ 
ple, in the case of multidimensional Einstein-Yang-Mills models with d extra dimensions, 
discussed in the next section, w = {d — l)/d |22j. Having in mind this class of models 
we take w ~ 0.5 -P 1. The exact value of this parameter does not affect our hnal result 
in a signihcant way. In the case of models obtained by dimensional reduction from mul¬ 
tidimensional theories with six extra dimensions, which can be motivated by the string 
theory, from Eq. m one gets Qm ~ 1-2. 

Fitting the Hubble law to the supernova dataset with a quadratic polynomial for G{z) 
— i.e., with only one free parameter, p 2 — does not give satisfactory results. To get a 


5 






log(z) 

Figure 1: The observational Hubble diagram of distant supernovae (dots with their corre¬ 
sponding error bars), the best £t curve to it in the model with the variable gravitational 
constant G{z) as given by Eq. (IT^ (solid line). 


better phenomenological approximation we consider a cubic polynomial as a parametriza- 
tion of the function G{z) in Eq. (IT^ . Varying p 2 and ps we obtained that the best fit to 
the Hubble diagram of SNela is achieved when the values of the parameters p2 and ps in 
Eq. (fT2|l are p2 ~ 0.34 and ps = —0.17, respectively. 

Fig. □ shows the observational Hubble diagram of distance moduli for SNela based 
on the data of [3] (their gold sample). Overplotted is the best £t curve for the predicted 
distance modulus 


=b\og di^z) + 2b+ (14) 

calculated in our model with the variable gravitational constant G{z), Eq. (I12|l . For these 
calculations we used the present Hubble parameter Hq = 63 km/s/Mpc. As it can be 
seen from Fig. ^ the agreement between the theoretical £t and the observations is quite 
good. 

The likelihood for the parameters p2 and pa is determined from the y^-analysis with 


X^(P2,P3) = ^ 
i 


[hth(2^i) - hobs]^ 

7 2 


(15) 
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Figure 2: Joint confidence intervals for the coefficients p 2 and of the polynomial £t to 
G{z) as given in Eq. (IT^ . 


where pth{z) is given by Eq. (UH), are the observational data for the distance moduli 
and cr^ are the uncertainties in the individual distance moduli. Dehnition (ITH|l of is 
analogous to the one in Ref. jS] . We would like to add that the best fit cubic polynomial 
G{z) yields = 199; for comparision the value of x^ obtained for the cosmic concordance 
model, which is x^ = 178- 

Fig. El shows the joint conhdence intervals for the £t to the SNela observational data. 
The analysis was done in the region of the parameters P 2 ,P 3 such that {z) > 0 for all 
redshifts in the interval 0 < z < 1.7. The bottom boundary of this region is seen in Fig. El 
To be more precise, for values of these parameters below the boundary the function g{z) 
dehned by Eq. o becomes negative as z approaches to z = 1.7. The best £t values of 
P 2 , P 3 are in fact rather close to this boundary. 

The plot of the cubic polynomial for G{z) with the coefficients found above is shown 
in Fig. El The phenomenological expression for G{z) based on the SNela data suggests 
that the value of the gravitational constant was larger in the past. To be more precise, 
G{z) > Go in the interval from z 0.03 to z ~ 1.97. It can be seen from Fig. El 
that for ^ > 0.03 as one moves towards larger redshifts the function G{z) hrst grows, 
reaches its maximum Gmax = 1.19 Go at z = 1.32 and then steadily decreases. Of course, 
the phenomenological expression for G{z) obtained here is approximate and presumably 
makes sense as far as the last, cubic term in Eq. m is smaller then the previous quadratic 
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Figure 3: The function G{z) represented by cubic polynomial of Eq. (I12jl with the coeffi¬ 
cients determined from the best £t to the SNela Hubble diagram. 

term, that is roughly for z < 2. 

The form obtained here for the function G{z) should be constrasted with existing 
constraints and/or compared with results on the variation of G in various cosmological 
models. As an illustration, in the following section G{z) will be considered as a phe¬ 
nomenological input in a multidimensional model and a prediction for the variation of the 
hne structure constant will be derived and analyzed. 

4 Variation of G and a in multidimensional models 

The time variation of the fundamental constants within models with extra dimensions 
has been considered in a number of papers j2I], IS1]^|3H], just to mention a few of them. 
Additionally, in ra the correlated variations of the hne structure constant and of the 
gravitational constant were analyzed. In this paper it was shown that in the framework of 
certain multidimensional models there exists a robust relation between the time derivatives 
of a and G which quite generically can be written as 


where the factor (3 is model dependent and, in general, may depend on the scale (or size) 



R of the space of extra dimensions. For the case of constant (3 a similar relation between 
the derivatives with respect to the redshift holds: 


Integrating this equation one gets 

Aa 

(y. OiQ y Crg J 

where = q;(0 ) is the present day value of the hne structure constant. 

To be specihc, let us consider the case of multidimensional Einstein-Yang-Mills the¬ 
ories. In this case it can be shown that /3 = 1. Using the cubic polynomial of Eq. m, 
with the coefficients determined in the previous section, the redshift dependence of a 
can be obtained and, from it, the behaviour of Aa/a predicted for such theory can be 
derived. In particular, it can be seen that a.t z = 0.5 this ratio is Aa/a ~ 0.06. This 
theoretical prediction is at odds with the known constraints on the rate of variation of 
the hne structure constant. More specihcally, the latest analysis of a Keck/Hires sample 
of quasar absorption lines using the many multiplet method gives 


--4 


a(z) — tto f G(z)\ 


(16) 




-= (-0.54 ± 0.12) X 10"® (18) 

a 

for in the range 0.5 < z < 3 [231 • Note, moreover, that these authors obtained that 
the value of Aa is negative, whereas the results derived from the Hubble diagram of 
SNela predict Aa > 0 for the same range of redshifts in the cosmological theories under 
consideration. Moreover, as already mentioned in Sect. 1, the observational results leading 
to a non-zero rate of variation of a have been challenged recently [23 EH- These studies 
provide much tighter constraints on the rate of variation of the hne structure constant. In 
particular, the constraint Aa/a = (—0.06±0.06) x 10“® was obtained, which is consistent 
with zero variation. In summary, neither the sign of our theoretical prediction nor its order 
of magnitude coincide with the above bounds. In fact, a similar conclusion was formulated 
in [221- 

The constraints studied before correspond to a redshift z ~ 0.5. The best local 
[z ~ 0.0) bound on the time variation of a is that obtained from the Oklo natural nuclear 
reactor: 

Aa , , 7 

-= 0.15 ± 1.05 X 10“^ 

a 

at z zz 0.15 iBg. The obtained phenomenological formula yields Aa/a ~ 0.006. Actually, 
as one can see from Eq. (ini), for Aa/a to ht the Oklo constraint at z = 0.15 the value of 
the parameter [3 should be of the order of |/d| ~ 10“®, which is the same to say that a{z) 
must be practically independent of G{z). 


5 Conclusions and discussion 

We have studied the possibility of htting the observational Hubble diagram of SNela 
assuming cosmological models of a flat universe without cosmological constant but with 
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a varying gravitational “constant”, G{z). The function G{z) was represented by a cubic 
polynomial, Eq. m, parametrized with three coefficients, pi, p 2 and p^. The linear order 
coefficient was hxed by the constraints on the present day rate of the time variation of 
G. The other two coefficients were determined from the best £t to the Hubble diagram. 
Finally we arrived at the following phenomenological expression: 

G{z) = Go (l - O.Olz + 0.34^2 - 0.17z^) , (19) 

where Go = G(0) is the present day value of the gravitational constant. It is important to 
mention that the sharp boundary at the bottom of the joint conhdence intervals in Fig. |21 
is due to the restriction that no negative values of H‘^{z) occur in the range of redshifts 
0 < 2; < 1.7. Values of ps smaller that those of the bottom boundary produce g{z) <0 
(see Eq. dH])) for some values of 2; in this interval and are consequently discarded. 

In this paper we limited ourselves to a cubic polynomial for G{z). Such choice of 
approximation is motivated by the results on two-parametric descriptions of astrophysical 
characteristics obtained from the same datasets in i> i and [21], which suggest that 
the available data does not allow a good determination of higher order parameters. The 
phenomenological determination of G(2;) from SNela data suggests that the value of the 
gravitational constant was higher than the present day one for 0.03 < 2; < 1.97 with the 
maximal value Gmax = 1.19 Go reached at 2; = 1.32. This change from the growing to 
the decreasing behaviour of G(2;) is a manifestation of the change in the behaviour of the 
observational data with the redshift. The latter feature was studied in Ref. |2j in terms 
of the deceleration q[z) represented by a linear polynomial with two parameters. Using 
the observational data for distant SNela the change of sign of q{z) at z = 0.46 ± 0.13 
was discovered. This is interpreted as an indication of the change from the epoch of 
acceleration of the evolution of the Universe to the epoch of deceleration as 2; increases. 

The features of the approximation G{z) obtained here should dehnitely be compared, 
within the domain of its validity, with other cosmological and astrophysical bounds and 
restrictions. The discrepancy between the theoretical prediction for the variation of the 
hne structure constant obtained in the models with extra dimensions and the existing 
observational constraints, provided that the latter are solid and confirmed, indicates that 
either the multidimensional models considered here are phenomenologically unsatisfactory 
or the very hypothesis of the variability of G is not correct. 
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